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Abstract. We present the first UBV and uvby photomet- 
ric observations for the short period variable star HV Ur- 
sae Majoris classified as a field RRc variable. The observed 
differences between the consecutive minima and the lack 
of colour variations disagree with the RRc-classification 
and suggest the possible binary nature of HV UMa. In 
order to reveal the real physical status of this star, we 
took medium resolution (A/ A A w 11000) spectra in the 
red spectral region centered at 6600 A. Spectra obtained 
around the assumed quadratures clearly showed the pres- 
ence of the secondary component. 

An improved ephemeris calculated using our and Hip- 
parcos epoch photometry is Hel. JD m ; n = 2451346.743 ± 
0.001, P= <i 7107523(3). A radial velocity curve was de- 
termined by modelling the cores of Ha profiles with two 
Gaussian components. This approximative approach gave 
a spectroscopic mass ratio of q sp =0.19±0.03. A modified 
Lucy model containing a temperature excess of the sec- 
ondary was fitted to the V light curve. The obtained set of 
physical parameters together with the parallax measure- 
ment indicate that this binary lies far from the galactic 
plane, and the primary component is an evolved object, 
probably a subgiant or giant star. The large temperature 
excess of the secondary may suggest a poor thermal con- 
tact between the components due to a relatively recent 
formation of this contact system. 

Key words: stars: binaries: eclipsing - stars: fundamental 
parameters - stars: individual: HV UMa 



1. Introduction 

The first note on the possible light variability of HV Ur- 
sae Majoris (= HD 103576 = HIP 58157, (V) = 8.69, 

Send offprint requests to: l.kiss@physx.u-szeged.hu 

* Based on the data obtained at the David Dunlap Observa- 
tory, University of Toronto 
** Guest Observer, Sierra Nevada Observatory 



AV = 0.28, P = d 355385, 7r Hi pp = 3.12 ± 1.23 mas, 
^Hipp = 320^gQ pc) was published by Penston (1973) who 
gave an 'uncertain' mark to the range of V-magnitude 
("var? V=8. 60-8.83"). The periodic nature of the light 
variation was discovered by the Hipparcos satellite (ESA 
1997) and the star was classified as an RRc variable. There 
is a note in ESA (1997) about the possibility of a dou- 
ble period but no firm conclusion was drawn. ESA (1997) 
gives a spectral type A3, while Slettebak & Stock (1959) 
published A7. 

We started a long-term observational project of 
Stromgren photometry and spectroscopy of the newly dis- 
covered bright Hipparcos variables. The first results have 
already appeared in Kiss et al. (1999a, b). Since the pe- 
riod, spectral type and light curve do not exclude the pos- 
sibility of wrong classification, accurate determination of 
the fundamental physical parameters is highly desirable. 

The main aim of this paper is to present the first U BV 
and uvby photometry for HV UMa. Also, our radial veloc- 
ity measurements are the first time-resolved spectroscopic 
observations of this star to date. The paper is organised 
as follows: the observations are described in Sect. 2, Sect. 
3 deals with data analysis and the obtained physical pa- 
rameters, while a discussion of the results is given in Sect. 
4. A final list of conclusions is presented in Sect. 5. 

2. Observations 

2.1. Photometry 

The Stromgren uvby photometric observations were car- 
ried out on 9 nights in June, 1999, using the 0.9 m tele- 
scope at Sierra Nevada Observatory (Spain) equipped 
with a six-channel (uvby+/3) spectrograph photometer 
(Nielsen 1983). Earlier, UBV measurements were ob- 
tained on one single night in March, 1999, using the 0.4 m 
Cassegrain-type telescope of Szeged Observatory equipped 
with a single-channel Optec SSP-5A photometer. These 
observations covered only 5 hours and revealed a 0.2 mag 
variation between two consecutive minima. We carried 
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Table 1. The journal of observations 



Julian Date 


type 


2451263 


UBV 


2451309 


spectr. 


2451310 


spectr. 


2451340 


uvby 


2451341 


uvby 


2451342 


uvby 


2451343 


uvby 


2451345 


uvby 


2451346 


uvby 


2451347 


uvby 


2451349 


uvby 


2451350 


uvby 



0.10 
0.12 
0.14 
£~ 0.16 
0.18 
0.20 



HV UMa - Sierra Nevada, 1999 June 
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r 
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0.16 
0.18 
0.20 
0.22 
0.24 
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Fig. 1. The light and colour curves of HV UMa phased 
with the adopted cphemeris (see text) 



out differential photometry with respect to HD 103150 
(V=8.45, B - V=0.54, b - y=0.335, mi=0.149, ci=0.381 
mag). The overall accuracy of the standard transforma- 
tion is about ±0.01 mag for V, b — y and mi and ±0.02 
mag for c\. The light and colour curves were phased us- 
ing the corrected ephemeris (see below) and are plotted in 
Fig. 1. 

2.2. Spectroscopy 

The spectroscopic observations were carried out at David 
Dunlap Observatory with the Cassegrain spectrograph at- 




6520 6530 6540 6550 6560 6570 6580 6590 6600 

MA) 

Fig. 2. Sample spectra around the quadratures and min- 
ima. The presence of the secondary component is obvious. 



tached to the 74" telescope on two nights in May, 1999. 
The detector and the spectrograph setup were the same as 
used by Vinko et al. (1998). The resolving power (A/AA) 
was 11,000 and the signal-to-noise ratio reached 30-50, de- 
pending on the weather conditions. The spectra were cen- 
tered on 6600 A and reduced with standard IRAF tasks, 
including bias removal, flat-fielding, cosmic ray elimina- 
tion, aperture extraction (with the task doslit) and wave- 
length calibration. For the latter, two FeAr spectral lamp 
exposures were used, which were obtained before and af- 
ter every three stellar exposures. The sequence of obser- 
vations FeAr-var-var-var-FeAr was chosen because of the 
short period of HV UMa. Careful linear interpolation be- 
tween the two comparison spectra was applied in order 
to take into account the sub-pixel shifts of the three stel- 
lar spectra caused by the movement of the telescope. We 
chose an exposure time of 10 minutes, which corresponds 
to 0.01 in binary orbital phase, avoiding phase smearing 
of the radial velocity curve. The spectra were normalized 
to the continuum by fitting a cubic spline, omitting the 
region of Ha. 

Besides a few telluric features, only the Ha line could 
be detected with acceptable S/N ratio in our 200 A-wide 
spectra. At the phases of maximum light the Ha line 
exhibited significant broadening and an excess bump ap- 
peared on the wings alternating between the blue and the 
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Fig. 3. Hipparcos epoch photometry data phased with the 
finally adopted ephemeris (see text) 



red side (see Fig. 2). It can be interpreted most easily as 
the effect of a close companion star, therefore HV UMa is 
most probably a spectroscopic binary. Fig. 2 shows a few 
sample spectra with the calculated orbital phase indicated 
on the right side of every spectrum. 

3. Physical parameters 

3.1. Epoch and period 

The Hipparcos data suggested that the light variation of 
HV UMa can be described with a single period of 01*355385 
(ESA 1997). We observed only one moment of minimum 
(Hcl. JD = 2451346.388), but the consecutive minimum 
appeared to be slightly fainter, therefore, we adopted a 
doubled Hipparcos period as a first approach (0 d 71077) 
and shifted the observed time of minimum with d 355385 
to obtain the final epoch Hcl. JD min = 2451346.743 ± 
0.001. 

The next step was to refine the period. This was done 
by phasing Hipparcos epoch photometry with the newly 
determined epoch and the doubled Hipparcos-period. The 
resulting phase diagram showed a shift of A(f> « 0.1 
(=0 d 071). That shift was eliminated by recalculating the 
period until correct phase diagrams for both our and Hip- 
parcos data (Fig. 3) were reached. The resultant period is 
P = d 7107523(3). The fact that earlier Hipparcos data 
agree very well with our data suggests a quite stable pe- 
riod of HV UMa. 

3.2. Classification 

The shape of the light curve, i.e. the continuous light vari- 
ation and the very deep secondary minimum (almost as 
deep as the primary one) , the absence of significant colour 
variation, the appearance of the secondary line in the spec- 
tra at the quadrature phases all suggest that HV UMa is 
probably an eclipsing contact binary. This is confirmed by 
its low mass ratio and a consistent model of the light curve 
(see below). 

Comparing the light curve phased with the final epoch 
and period (Fig. 1) with the line profiles observed at 



phases of maximum light (Fig. 2) it is visible that the 
secondary line appears on the blue side at the <j> = 0.27 
quadrature phase that follows the deeper minimum, while 
this bump is rcdshiftcd at (j> = 0.77. This indicates that the 
smaller companion star approaches us after the primary 
minimum, therefore that minimum is due to an occulta- 
tion eclipse. The weak point of this analysis is that the 
light curve is not very well covered around the minima ci- 
ther by the Hipparcos data or our observations. The data 
obtained at Sierra Nevada have better inner precision (less 
scatter) than those provided by Hipparcos, and these data 
indicate that the minimum at (j> = 0.0 is slightly deeper. 
Therefore, we adopted this eclipse as primary minimum, 
but this needs further confirmation. If the deeper minim- 
mum is really due to an occultation eclipse then HV UMa 
is a so-called W-type contact binary. 

Contact binaries can contain early (O-B) or late (G- 
K) spectral type stars. The latter group is referred to as 
the W UMa stars, while the former is known as early-type 
contact systems. The colours of HV UMa indicate early 
F spectral type, therefore HV UMa is an "intermediate" 
type contact binary between the W UMa stars and the 
OB-type contact systems. The surface temperature and 
the line profile of the HV UMa system makes it similar to 
the known contact systems UZ Leo and CV Cyg (Vinko 
et al., 1996). 

3.3. Radial velocities and spectroscopic mass ratio 

Since the secondary component is only partly resolved, 
the radial velocities must be determined by modelling the 
individual line profiles in order to avoid blending effects 
(e.g. systematic decrease of the velocity amplitude). For 
this purpose we chose those spectra that were obtained 
around the quadratures. These show the presence of the 
secondary most clearly. One spectrum around light mini- 
mum was also modelled to test the applied method. 

Because the Ha profile is strongly affected by the Stark 
broadening and shows wide non-gaussian wings, we nor- 
malized the profiles to the surrounding continuum, and 
selected the lower part of the profiles below the 0.9 in- 
tensity value. We fitted two individual Gaussian profiles 
to the line cores adjusting the amplitudes, FWHM values 
and line core positions. The initial values of these param- 
eters were estimated from two spectra very close to the 
quadratures (0=0.25 and 0.77). The FWHM converged 
very quickly to the final values, being 7.6 A and 4.0 A for 
the primary and secondary components, respectively. Line 
depths changed slightly from spectrum to spectrum, as the 
contributions are phase-dependent, resulting in 0.29-0.30 
for the primary and 0.07-0.10 for the secondary (note, 
that these values mean line depths below 0.9 normalized 
intensity). The fitted line core positions resulted in the ra- 
dial velocity variations for both components. Sample spec- 
tra with the fitted profile are shown in Fig. 4, while the 
radial velocities are presented in Table 2. The estimated 
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accuracy of the individual velocities is about 5 km s _1 for 
the primary, and 10 km s^ 1 for the secondary, which is 
mainly determined by the resolution of the line core in 
wavelength. The velocity amplitudes resulting from this 
method are 47±1.5 km s _1 and 254±10 km s _1 , where 
the uncertainties are due to the random errors caused by 
the observational scatter. The corresponding mass-ratio is 
<7s P = m soc /m pri = 0.185 ± 0.01. 

However, as was also pointed out by the referee, this 
kind of velocity determination may contain a large amount 
of systematic error, mainly due to the assumed Gaussian 
shape of the individual line profiles. The intrinsic Ha pro- 
files of the components of HV UMa are probably quite 
different from Gaussian, therefore this approach can be 
considered as only the first approximation for extracting 
the radial velocities from the Ha profiles. The major part 
of the systematic error is governed by the shape of the 
wing of the primary component's model profile on the side 
where the secondary star appears (blueward at <f> = 0.25 
and redward at </> = 0.75 phases). It is well visible in Fig. 
4 that the position of the secondary line is shifted to- 
ward larger velocities with respect to the position of the 
"hump" on the observed profile, due to the increasing con- 
tribution of the primary line toward the main minimum 
of the combined line. If the primary line was steeper on 
the side where the secondary line exists, overlapping the 
secondary by a smaller amount, then the secondary line 
would be less shifted, thus, its position would be closer 
to the local hump on the observed profile, resulting in 
a smaller radial velocity of the secondary. On the other 
hand, a shallower secondary profile would give us system- 
atically higher velocities due to the same reason. 

In order to estimate the amount of this kind of sys- 
tematic error, we simply determined the positions of the 
two local minima (the main minimum and the secondary's 
hump) on the profiles observed around quadratures (four 
spectra around <p — 0.25 and two around <fi — 0.75) when 
the presence of the hump appeared to be most promi- 
nent. This was done interactively, by eye, plotting the line 
profiles on the computer screen, which again introduced 
some subjectivity into the procedure, but it is stressed 
that this is done only for estimating the errors of the veloc- 
ities and not for obtaining their actual values. Of course, 
the velocities of the secondary measured in this way were 
systematically smaller than those obtained by the Gaus- 
sian fitting. The velocities of the primary were almost the 
same, as could be expected. The total amplitude turned 
out to be K' = 280 km s _1 , while the mass ratio changed 
to q' = 0.22. Comparing these values with the results of 
the Gaussian fitting, we conclude that the errors of the 
radial velocity amplitude and the spectroscopic mass ra- 
tio (both random and systematic) are approximately ±23 
km s -1 and ±0.03, respectively. The finally adopted pa- 
rameters determined spectroscopically, together with their 
errors are collected in Table 3. It is important to note that 
the mass ratio can be refined by modelling the light curve 



Table 2. The observed heliocentric radial velocities ob- 
tained by the Gaussian fit. The velocity resolution is about 
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Hel. JD 
2400000+ 





V Tad {prim.) Kad (sec.) 
[km/s] [km/s] 


51309.6206 
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-47 


251 


51309.6396 





.80 


-51 


242 


51309.6471 
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-42 


246 


51310.6411 





.21 


41 


-252 


51310.6485 





.22 


45 


-243 


51310.6562 





.23 


50 


-243 


51310.6643 





.24 


54 
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51310.6719 





.25 


50 


-243 


51310.6794 





.26 


50 


-243 


51310.6893 





.27 


45 


-238 


51310.6964 





.28 


41 


-252 


51310.7037 





.29 


45 


-247 


51310.7436 





.35 


41 


-233 


51310.8423 





.49 
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(Sect. 3.4), but the total velocity amplitude is tied only to 
the spectroscopic data, thus, its uncertainty will directly 
appear in the absolute parameters of the system. 

3.4- Light curve modelling 

The V-light curve was synthesized with the computer code 
BINSYN described briefly in Vinko ct al. (1996). This code 
is based on the usual Roche-model characterized by the ge- 
ometric parameters <;ph (photometric mass-ratio), F (fill- 
out) and i (orbital inclination). The relative depth of the 
eclipses were modelled introducing the relative tempera- 
ture excess of the secondary X — (T scc — T pr i)/T pr i (hot- 
secondary model). Because the primary minimum turned 
out to be due to occultation, the phases were shifted by 
0.5 assigning 4> = 0.0 to the transit eclipse (built-in default 
in BINSYN). 

First, the effective temperature of the primary com- 
ponent was estimated based on synthetic colour grids by 
Kurucz (1993) and the observed mean Stromgrcn colour 
indices ((b — y) = 0.19 mag, (mi) = 0.15 mag, (ci) = 0.77 
mag), resulting in T off = 7300±200 K and log .9 = 4.0±0.3 
dex (assuming E(B — V) = and solar chemical abun- 
dance). The interstellar reddening in the direction and 
distance of HV UMa is expected to be small, because this 
variable lies far from the galactic plane. TU UMa, an RR 
Lyr variable lying 17° SE from HV UMa also has a negli- 
gible colour excess (Liu & Janes, 1989). 

Other parameters necessary for modelling the binary 
star were as follows. A linear limb-darkening law with co- 
efficient u — 0.61 was adopted from tables of Al-Naimiy 
(1977). The gravity darkening exponent and the bolomet- 
ric albedo were chosen at their usual values for radiative 
atmospheres: (3 — 0.25 and A = 1.0. All these parameters 
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Fig. 4. The observed (crosses) and fitted (solid lines) 
Gaussian line profiles at two quadrature phases. 



were kept fixed during the solution for the best light-curve 
model. 

The light-curve fitting was computed using a con- 
trolled random search method, the so-called Price algo- 
rithm (Barone et al., 1990, Vinko et al., 1996). The op- 
timized parameters were g p h, F, i and X. The best solu- 
tion was searched for in the following parameter domains: 
0.05 < q < 0.5, 1.01 < F < 1.99, 50 < i < 90 and 
< -0.2 < X < +0.2. The fit quickly converged to low 
inclination and low mass ratio values that were expected 
from the shallow eclipses (AV « 0.3 mag) and the small 
spectroscopic mass ratio (q sp — 0.19). Also, it turned out 
that there are strong correlations between the optimized 
parameters. Due to this correlation, the physical param- 
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Fig. 5. The observed (crosses) and calculated Ha profiles 
at 0=0.25. The spectrum was calculated by an ATLAS9 
code adopting T e g — 7000 K, log 5=4.0 and full radial ve- 
locity amplitude of 300 km s _1 . The bottom panel shows 
the residual intensities for the adopted fit compared with 
two other models (6750 K and 7250 K). The overall agree- 
ment is the best for T c tf = 7000 K. The sharp absorption 
features are atmospheric telluric lines. 



ctcrs determined from the light curve fitting cannot be 
considered as a unique solution: certain parameter com- 
binations describe the light curve almost equally well. In 
these cases the combination of photometric and spectro- 
scopic information is very important: one can use the spec- 
troscopic data to determine a consistent set of physical pa- 
rameters that gives the best model fitting all the available 
data. 

In order to combine the photometric and spectroscopic 
information and find a consistent description of the sys- 
tem, we modelled the observed Ha line profiles using the 
parameters from the light curve fitting. The model line 
profiles were computed by convolving an intrinsic Ha pro- 
file of a non-rotating star with the Doppler-broadening 
profile of the contact binary. The broadening profiles were 
calculated with the WUMA4 code (Rucinski, 1973). For 
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Fig. 6. Fitted light curve of HV UMa. Note a 0.5 shift in 
phase to get a transit eclipse at </>=0. 



the determination of the intrinsic Ha profile we used Ku- 
rucz's ATLAS9 code modified by John Lester. This ap- 
proach, however, has some limitations, because the Ha line 
is strongly affected by NLTE-mechanisms, therefore the 
ATLAS9-model profile will be somewhat different from the 
real intrinsic profile, especially near the line core. Thus, 
only a crude comparison between the modelled and ob- 
served line profiles was possible, neglecting the differences 
in the line core. 

However, it turned out that the originally adopted ef- 
fective temperature T c s = 7300 K was too high. With this 
temperature, the Stark-wings of the Ha line are so strong 
that they overwhelm the rotational broadening, produc- 
ing much wider line profiles than observed. Thus, we re- 
duced the value of the effective temperature until satisfac- 
tory agreement was found between the observations and 
the broadened model profiles. The resulting fit is plotted 
in Fig. 5 together with the observed line at <j> = 0.25. 
The residuals of three model profiles are also shown. It 
can be seen that at T e g = 6750 K the wings can be fit- 
ted quite well, but the computed line core is too shallow. 
On the other hand, if T e g — 7250 K is used, the com- 
puted line core agrees better, but the wings are wider. 
Thus, T c g = 7000 K was accepted as a compromise. The 
Stromgren-colours of the Kurucz-grid for this temperature 
are almost the same as for the originally adopted 7300 K, 
therefore this temperature is still consistent with the pho- 
tometric colours. 

The light curve modelling was recomputed with T e g = 
7000 K. The optimized parameters changed only slightly, 
resulted in a slightly smaller mass ratio and a slightly 
higher temperature excess. Table 3 shows the final set of 
physical parameters, while the fitted light curve is visible 
in Fig. 6. The distribution of the random points around the 
X 2 -minimum in the parameter space is plotted in Fig. 7. 




I_i I i I i I I 1 75 L_, I i I i I I 

0.182 0.184 0.186 57.2 57.3 57.4 

q i 

Fig. 7. Distribution of random points around the mini- 
mum of x 2 ■ 



The structure of the sub-spaces in this diagram indicates 
the correlation between the different parameters. 

Another light curve model was computed in order to 
test the effect of the gravity darkening and reflection pa- 
rameters that were originally fixed as if the atmosphere 
was radiative. The model with their "convective" values 
(3 = 0.08 and A = 0.5 resulted in an even larger tem- 
perature excess than in the radiative case. Because the 
temperature excess of the secondary is only a "correc- 
tion" parameter in the light curve solution and it may 
not mean real temperature difference, it would be difficult 
to explain physically a very large value of the temperature 
excess that still does not cause significant colour variation. 
Thus, we adopted the model with radiative atmospheric 
parameters as our final solution, and this model is listed in 
Table 3. Note that the errors of the fitted parameters (3rd 
column) are difficult to estimate, because of the parameter 
correlation. We monitored the behaviour of the \ 2 func- 
tion during the optimization and assigned uncertainties 
to each parameter according to the spread of the random 
points for which \ 2 < 0.5. This criterion defines those so- 
lutions when the fitted curve runs well within the error bar 
at each measured normal point, giving a feasible fit to all 
observations. The parameter correlation means that the 
uncertainties are also not independent of each other: e.g. 
slightly decreasing q forces increasing fori (see Fig. 7). 
The uncertainties of the calculated parameters (3rd panel 
in Table 3) were estimated assuming a ±23 km s" 1 error 
in the radial velocity amplitude K. 

Because of the correlation between the optimized pa- 
rameters, it is very important to check whether the radial 
velocities calculated from the model match the observed 
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Table 3. Physical parameters of HV UMa. 



Fitted parameters 


Value 


o 


spectroscopy 


K (km s" 1 ) 


300 


23 


Qsp 


0.19 


0.03 


Teff (K) 


7000 


200 


photometry 




0.184 


0.05 


F 


1.77 


0.15 


i 


57.3 


0.4 


X 


0.13 


0.03 


Calculated parameters 


Value 


G 


a (10° km) 


3.48 


0.25 


Mi(M ) 


2.8 


0.6 


M 2 (M ) 


0.5 


0.17 


Ri(Re) 


2.62 


0.25 


i?2(-R©) 


1.18 


0.16 


ppri (g cm -3 ) 


0.2 


0.05 
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Fig. 8. The radial velocity curves of the primary and sec- 
ondary components. Symbols correspond to the directly 
measured values, while solid lines denote the calculated 
ones from the photometric model. 



velocities. This comparison is plotted in Fig. 8 where an 
almost perfect agreement can be seen. The low mass ratio 
results in the distortion of the sinusoidal velocity curves, 
which is a well-known effect in close binaries. 

4. Discussion 

The spectroscopic detection of the secondary component 
and the success of modelling the light- and velocity curves 
as well as the Ha line profile supports our conclusion that 
HV UMa is not a short-period RR Lyr variable, but an 
eclipsing binary system. The physical parameters listed in 
Table 3 give a deep contact configuration of this binary, 
explaining naturally the lack of significant colour change 
during the light variation cycle, which would be peculiar in 



a pulsating variable. Since the primary minimum is due to 
occultation, HV UMa is a so-called W-type contact binary 
(systems exhibiting transit eclipses as primary minima are 
called A- type). 

It is well known that contact binaries can be formed ei- 
ther from hot, early- type stars or cool, late- type stars, the 
latter representing the class of W UMa-type variables (see 
e.g. Rucinski, 1993; Figueiredo et al., 1994 for reviews). 
The surface temperatures of W UMa-systems are usually 
below 7000 K, while the temperatures of early- type con- 
tact binaries show a much wider range, between 10,000 
and 40,000 K. Thus, HV UMa falls into the temperature 
regime between 7000 and 10,000 K that is relatively rarely 
occupied by contact binaries. Systems like HV UMa may 
represent a transition population between early-type and 
late-type contact binaries (although even the existence of 
such transition is questionable). These systems lie on the 
boundary between radiative and convective atmospheres, 
at about late A - early F spectral types. Detailed studies of 
such systems would be important, because the formation 
and the structure of early-type systems having radiative 
envelopes and late-type contact binaries having convective 
envelopes is probably quite different. 

The period value P w 0.711 day also indicates that 
HV UMa is not a typical contact system. Recent statisti- 
cal studies based on the data from the OGLE microlens- 
ing survey (Rucinski, 1998 and references therein) have 
shown that the number of contact systems strongly de- 
creases above P = 0.7 day, and there is a well-defined limit 
at P w 1.5 days. Therefore, HV UMa is a member of the 
relatively rare "longer-period" contact binaries, although 
a few contact systems with P > 2 days definitely exist, at 
least close to the galactic bulge. On the other hand, there 
exists a period-colour relation of "normal" W UMa-stars 
with period P < 1 day, which indicates that longer period 
systems have bluer colour. On the plot of the empirical 
logP — (b — y)o relation (Rucinski, 1983) the position of 
HV UMa, using the mean colour (b — y) = (b — y)o = 0.19 
(Fig.l, assuming zero reddening) and our revised period 
logP = —0.1483 (Sect. 3.1), is close to the upper bound- 
ary of the relation, suggesting an atypical, but not peculiar 
contact system. 

A more recent logP — (B — V)o diagram based 
on Hipparcos-parallaxes has been published by Rucinski 
(1997, see his Fig.2). The position of HV UMa on this 
diagram was calculated assuming E(B — V) = and 
(b - y)/(B - V) = 0.7, resulting in (B - V) « 0.27. 
These data show that the position of HV UMa is entirely 
consistent with that of involving the older logP — (b — y) 
relation, being close to the blue short-period envelope 
(BSPE, Rucinski, 1997), but the system is definitely red- 
der than the upper limit defined by the BSPE, consis- 
tently with other contact binaries. Also, a rough com- 
parison of HV UMa with other contact binaries on the 
logP — (V — J)o diagram based on OGLP-photometry 
(Rucinski, 1998) strengthens the status of HV UMa out- 
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lined above, again, being closer to the BSPE than other 
systems with similar period, although the lack of observed 
(V — I) colour of HV UMa limits the reliability of this 
comparison at present. It can be concluded that all avail- 
able measurements and pieces of information consistently 
support the contact binary nature of HV UMa. 

The new physical parameters collected in Table 3, to- 
gether with the Hipparcos-parallax (ir = 3.12 ± 1.23 mas) 
enable us to estimate the evolutionary status of HV UMa, 
provided it is indeed a contact binary with those pa- 
rameters. The coordinates and the distance based on the 
Hipparcos-parallax indicate that HV UMa is a halo ob- 
ject: its distance from the galactic plane is z « 300 pc, 
which means that [Fe/H] = (assumed during the anal- 
ysis of the line profiles and the colour indices) may not 
be true. Because the spectroscopic observations presented 
in this paper have limited spectral range (AA = 200 A 
around 6600 A), and this region does not contain signif- 
icant metallic lines in early-type stars, the spectroscopic 
derivation of [Fe/H] was not possible. Thus, because of the 
lack of further information we assumed [Fe/H] = 0, which 
is not impossible for halo objects, but [Fe/H] < is also 
likely. 

The referee suggested the possibility that the low metal 
content of HV UMa critically affects the Stark broadening 
of the hydrogen lines, thus, significantly influencing the 
temperature derived from the Ha profile in Sect. 3.4. We 
investigated this effect in detail using the pre-computcd 
Ha profiles by Kurucz (1979) including Stark broaden- 
ing. In the left panel of Fig. 9 the dependence of the Ha 
equivalent width on metallicity is plotted. For each tem- 
perature, two model sequences for \ogg = 4.0 and 3.5 are 
shown. It can be seen that the decrease of the metallicity 
indeed affects the strength of Ha, but in this temperature 
range the variation of the equivalent width is governed 
mainly by the change of the effective temperature. The 
gravity (pressure) dependence is very weak. Because the 
equivalent width of the broad Ho: line strongly depends 
on the strength of the Stark wings, it is expected that the 
wings of the Ha profile presented in Fig. 5 (corresponding 
to [A/H] = 0) arc not affected very largely by the possible 
lower metallicity of HV UMa, thus, the derived temper- 
ature T c ff = 7000 K is only slightly dependent on metal 
content. This is illustrated in the right panel of Fig. 9, 
where two model profiles corresponding to two different 
temperatures (T cff = 7000 and 7500 K) and [A/H] = -1.0 
i.e. significantly lower metallicity than assumed in the pre- 
vious section are presented together with the observed line 
profile at quadrature. It can be seen that the T e s = 7500 
K model still gives a broader line profile than observed, 
while the T e g = 7000 K model results in a much better 
agreement in the wings, very similar to the case of solar 
metallicity presented in Fig. 5. Note, however, that the 
lower metal content causes a less deep line core of Ha, 
thus, the problem of fitting the whole Ha line is exagger- 
ated when the effect of metallicity is taken into account. 



Nevertheless, it is concluded that the T cS = 7000 ± 200 
km s" 1 temperature derived from the wings of Ha proba- 
bly does not contain a significant systematic error due to 
the unknown metallicity of HV UMa. 

As was mentioned in Sect. 3.4, the colour excess of 
HV UMa is E(B - V) w 0. This is supported by the ef- 
fective temperature derived spectroscopically (discussed 
above) and photometrically (from observed and tabulated 
Stromgren indices), because both methods resulted in a 
consistent value. The negligible reddening is also in agree- 
ment with the statement that HV UMa belongs to the 
halo population. 

At first glance, the absolute geometric parameters col- 
lected in Table 3 would indicate that the HV UMa system 
consists of main-sequence components: both stars have 
log g — 4.0 and the mass and radius values of the pri- 
mary are also similar to those of a main sequence star (the 
secondary is oversized in relation to its mass, typical of 
contact systems). However, the surface temperatures and 
luminosities indicate that HV UMa is probably an evolved 
object. First, the combined absolute magnitude of the sys- 
tem based on parallax measurement and E(B — V) = 
results in My — 1.0 ± 0.8 mag, where the large error is 
due to the uncertainty of the Hipparcos-parallax. Using 
tabulated bolometric corrections, the total luminosity of 
the system is Lt = 30 ± 20L Q . Second, the luminosities 
of the components are L; = in Rf oT 4 ff = 15.4L Q and 
2.9-Lq for the primary and secondary, respectively, giving 
Lt = 18.3L0 for the combined luminosity, which is within 
the error of the distance-based total luminosity estimated 
above. However, both of these luminosities are much less 
than the expected luminosity L w 60 ± IOLq of a main 
sequence star with M « 3M© (Lang, 1991). Moreover, 
this kind of main sequence star would have T ff = 9500 K, 
much higher than the surface temperature of HV UMa. 
Therefore, the primary component of HV UMa is too cool 
and too faint for its mass if it is assumed to be a main 
sequence object. 

The agreement with a class III giant star having L « 
30 ± IOLq for M w 3M Q is much better. The tempera- 
ture of such giant star is T e g w 7400 ± 300 K which is 
not very far from the surface temperature of HV UMa. 
Taking into account the energy transfer between the com- 
ponents in the contact binary (assuming that the total 
luminosity of the system is due to the energy production 
of only the more massive primary component), the cor- 
rected effective temperature of the primary component is 
7\,corr = 7300 ± 100 K. The radius and the surface gravity 
of this giant star, R = 2.9R@ and logg « 3.8, also agrees 
well with the derived parameters of HV UMa. Therefore, 
the comparison of empirical and theoretical values of the 
physical parameters suggests that the primary component 
of HV UMa is an evolved object, probably a IV-III class 
subgiant, or giant star. Because W UMa stars are gener- 
ally accepted to belong to the old disk population (e.g. 
Rucinski, 1993, 1998), it is reasonable that a long-period 
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Fig. 9. Effect of metallicity on the Ha line profile. Left panel: the dependence of the Ha equivalent width on the 
metallicity for three temperatures (labeled) and two \ogg values (4 and 3.5) for each temperature. Right panel: 
comparison of the Ha of low metallicity models (for T e g = 7000 and 7500 K) with the observed profile. 



contact system, containing a more massive primary than 
most of other W UMa systems, is significantly evolved 
from the main sequence. Therefore, the evolved status of 
HV UMa qualitatively agrees with the age of other contact 
binaries. 

It is interesting to compare the direct empirical ab- 
solute magnitude of HV UMa derived above (My = 1.0 
mag) with the prediction of the period-colour-luminosity 
relation of W UMa stars calibrated by Rucinski & Duer- 
beck (1997) as M v = 0.10 + 3.08(5 - V) - 4.42 log P. 
Using the same estimated (B — V)o index as above, 
the predicted absolute magnitude for HV UMa becomes 
My = 1.59 ± 0.35 mag, which agrees with the empirical 
value within the errors. Note, that the deviation of some 
of the calibrating W UMa stars in the sample of Rucinski 
& Ducrbeck (1997) from the value predicted by this rela- 
tion is as large as 0.5 — 0.7 mag (see Fig. 4 in Rucinski & 
Ducrbeck 1997), therefore the difference between the ob- 
served and the predicted absolute magnitude of HV UMa 
does not make this system discrepant with respect to other 
contact binaries. On the other hand, it is a bit surprising 
that the relation that is mainly based on main sequence 
objects gave such a good prediction for the more evolved 
HV UMa system. This agreement is probably limited to 
the particular range on the HR-diagram close to the po- 
sition of HV UMa, and may not hold on for more evolved 
systems with P > 1 day. Very few known contact systems 
exist above the P = 1 day period value, as recently dis- 
cussed by Rucinski (1998), this lack of systems also gives a 
natural limit for the applicability of this relation for longer 
periods. 

The separation of the components in the HV UMa sys- 
tem and the evolved physical state of the primary may 
suggest that this contact system formed during a case 
B mass transfer. This may also give a reasonable expla- 
nation for the poor thermal contact AT = 900 K be- 
tween the components. Model computations of the for- 
mation of contact binaries via evolution induced mass 



transfer from the more massive component (Sarna & Fe- 
dorova, 1989) predict large amount of temperature excess 
(AT « 2000 - 3000 K) at the moment of reaching the 
contact configuration. The result that the eclipse depths 
of HV UMa can be modelled with only such high temper- 
ature excess may indicate that this contact system formed 
only recently and did not have enough time to reach bet- 
ter thermal contact. Note, that the temperature excess in 
W UMa-type contact binaries is usually considered un- 
physical, because the lack of the colour index variation 
suggests very good thermal contact for late-type stars. 
The physically consistent model of the eclipse depths of 
W UMa-stars contains large starspots on the surface of 
one or both components (e.g. Hendry et al., 1992). How- 
ever, in the case of HV UMa with T eff >= 7000 K, the 
presence of such starspots is less likely, thus, the eclipse 
depths of this system may indeed mean a 900 K temper- 
ature difference between the secondary and the primary. 

5. Conclusions 

The new results presented in this paper can be summa- 
rized as follows. 

1. We reported the first uvby photometric and medium 
resolution spectroscopic observations of the short-period 
variable star HV Ursae Majoris. Contrary to the RRc clas- 
sification given by ESA (1997), the star turned out to be 
a new contact binary with early-type components. 

2. An improved ephemcris was determined us- 
ing our and Hipparcos epoch photometric data: 
Hel.JD min =2451346.743+0.7107523(3)-E. There is 
no indication of changing period over almost 10 years. 

3. A radial velocity curve was measured directly from the 
Ha profiles by two-component Gaussian fitting of the line 
core regions in spectra recorded around the quadratures. 
We calculated a spectroscopic mass ratio of 0.19±0.03. 

4. The effective temperature and the surface gravity were 
determined using the mean Stromgren indices, synthetic 
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colours of Kurucz (1993) and theoretical (ATLAS9) line 
profiles giving T ef f = 7000 ± 200 K and log g = 4.0 ± 0.3. 
The light curve modelling resulted in a complete set of 
physical parameters. 

5. The physical parameters of HV UMa together with the 
parallax measurement indicate that this binary is situated 
far from the galactic plane, and the primary component is 
an evolved object, probably a subgiant or giant star. The 
large temperature excess of the secondary may be indica- 
tive of a poor thermal contact between the components 
due to a relatively recent formation of this contact system 
via case B mass transfer. 
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